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We investigate the possibility of probing the CP-violating htt coupling in the process e"'"e“ —>■ /17 
at the future high luminosity e^e“ colliders. Our numerical results show that the cross section for 
this process can be significantly increased for the allowed CP phase ^ and center of mass energy. For 
example the cross section is about 10 times of that in the standard model (SM) for ^/s = 350 GeV 
and ^ = Stt/S (see text for ^ definition). The simulation for the signal process e^e“ —>■ ^,7 —>■ blry and 
its backgrounds shows that the signal significance can reach about 5cr and more than 2.1a for ^/s = 
350 GeV, 500 GeV respectively, with the integrated luminosity £ = 3 ab“^ and ^ £ [7r/2, 37r/5]. 
For C — IQ ab“^, the signal significance can be greater than 5 ct for y/s = 350 GeV and about 4 ct for 
7s = 500 GeV with the CP phase ^ G [7r/2,37r/5]. Besides the cross section enhancement, the CP- 
violating hit coupling will induce a forward-backward asymmetry Afb which is absent in the SM 
and is a clear signal of new CP violation. Compared with the Afb in the Higgs decay h —>■ l'^l~'y, 
the Afb can be greatly enhanced in the production process. For example Afb can reach -0.55 for 
^ = 7r/4 and 7s = 500 GeV. Due to the large backgrounds, the significance of the expected Afb 
can be only observed at 1.68 (t with C — IQ ab“^ and 7® = 500 GeV. It is essential to trigger the 
single photon in the final state to separate the bottom jets arising from scalar or vector bosons, in 
order to isolate the signal from the backgrounds more efficiently. 
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I. INTRODUCTION 


The observation of the 125 GeV Higgs boson at the Large Hadron Collider (LHC) (TJ H] marked 
a milestone in particle physics. Consequently detailed measurements of the discovered Higgs boson 
properties have become one of the main priorities of the LHC and future colliders, in order to verify 
whether it is the SM Higgs boson or not. In fact, there are various motivations for new physics beyond 
the SM (BSM) [3]. In the SM, charge conjugation-parity (CP) violation is described by a 3 x 3 quark 
mixing matrix, the well-known Cabibbo-Kobayashi-Maskawa (CKM) matrix [HE], with a single complex 
phase in the gauge sector, while the Higgs sector is CP-conserved. However, the CP violation in the SM 
can not account for the origin of the baryon asymmetry of the universe (BAU) [SJIZ] which is characterized 
by the baryon-to-entropy ratio ns/s ~ (8.59 ± 0.10) x 10“^^ [HIIS]- Thus it is necessary to look for new 
sources of CP violation. 

Phenomenologically there are already many works on studying the CP-violating Higgs couplings |10L 
[35] . In general, the CP-odd gauge-Higgs couplings are generated from higher dimension operators HO] 
and usually smaller compared to the CP-even couplings. However /177, /1Z7 couplings are not present 
at the tree-level, so CP-even and CP-odd components can be of the similar magnitude. Especially the 
/177, /1Z7 couplings can be induced by the CP-even and CP-odd Yukawa couplings Hi!, which are not 
suppressed and comparably large (33] [35] . Due to the interference of the fermion loops and W boson 
loops, the /177, hZj couplings are sensitive to the Higgs coupling to t quarks [numsa. 

The forward-backward asymmetry (Afb) as a consequence of parity violation has been studied exten¬ 
sively in Z physics at the LEP |33|- In [HHiiiiii, the authors proposed to measure the CP violation in 
Higgs decays using this observable. Since the intermediate Z boson is dominantly on-shell but the photon 
is virtual and the imaginary part is proportional to the ■j — Z interference term, the Afb oi the process 
h — 7 , Z+U 7 is suppressed by a factor oiTz/mz ^ 3% [HJUT]. However, the Afb will be greatly enhanced 
in the production process e“'"e“ —>■ /ly, where the /177 and hZj couplings can have large imaginary parts 
(due to resonance effects) |33|- In this paper, we will discuss the possibility of probing the CP property 
of hit coupling in e“''e“ —>■ /ly at the future e“''e“ colliders. 

The paper is organized as follows. In Sec.jHj we will present the one-loop amplitude of e“'"e“ —)■ h"f and 
the analytical formula of Afb as well as the CP symmetry of the helicity amplitude. In Sec. jnTj we give 
the numerical results with the help of FeynArts/FormCalc/LoopTools |3SJ[3S] for some benchmark center- 
of-mass (c.m.) energies and CP phases. In Sec. IV the collider simulation of signal and backgrounds 
is presented. The analytical expressions of scalar functions, W boson loop functions and box diagrams 
contributions and the derivation of the significance of Afb are collected in the appendices. 


II. FORMALISM 

We first parameterize the CP-violating hit coupling as |33| 

TTit — 

Chtt = - Ktt{cos^t + ism(tl5)t, ( 1 ) 

V 

where v = 246 GeV is the vacuum expectation value (vev) of the Higgs field, Kt G M. describes the 
magnitude of hit interaction and ^ S (—7r,7r]. The CP-even and CP-odd cases correspond to ^ = 0 
and ^ = 7r/2, respectively. In particular, the SM Higgs boson has Kt = I and ^ = 0. The signal strengths 
measured in gg ^ h and h ^ -fj have constrained the CP phase |^| < 37r/5 at 95% C.L. |23l|27]^. On 
the other hand, the electric dipole moment (EDM) of the electron de constrains the upper limit on the 
CP-odd coupling |Ktsin^| < 0.01 with the assumption of SM-like hee coupling [3n]. However, this can 
be evaded in various new physics (NP) models |3H||151I15]. We will suppose Kt = 1 and ^ £ [0, 37r/5] in 
this paper. 

The process e+e“ —>■ /ly has been investigated in the SM [5UH53] and beyond the SM |311[53H55] . The 
tree-level amplitude is suppressed by the small hee Yukawa coupling and can be safely ignored EB. Fig.|I| 


^ ATLAS Collaboration has performed a global analysis of the Higgs couplings HZ], which may give a stronger constraint 
on the htt coupling, however it depends on the assumption of the hgg couplings as well so the conclusions in our paper 
will not change. 




3 


shows the representative diagrams for e+e” —>■ h'y at one-loop level^, where only the htt Yukawa coupling 
is assumed to be CP-violated and contained in /lyy and hZ^ vertex diagrams and j/Z — h mixing 
diagrams. The latter is the higher order correction to hee coupling, which is also neglected. The /lyy 
and hZ^ effective couplings take the form as |53] 

( 2 ) 

where q is the four-momentum of the intermediate propagator V = 'y, Z and the coefficients are 

generally complex, 


G] 

Gf 


^ e mw ^pi,w _ 

sw 







(3) 


where y = 1,6 for i = 1 while y = — 1 for z = 3. Qf, rrif and Nc are the electric charge, mass of the 
fermion and color factor, respectively, and the vector part of Zff coupling gy = — 2 sin 9wQf with 

the third component weak isospin the left-handed fermion /. The Higgs boson is C-even 

and the photon is C-odd, only the C-odd (vector) part of the Zff coupling gy contributes [301111] even 
though the htt coupling is CP-violated [ITj. The fermion-loop contribution functions fI , f( and fI are 
given by (^/ = 0 for / 7 ^ t) 

fI = -—cos - s){-4m}+ml-s)Go + 2s{Bois)-Boiml)) + 2{s-ml)], 

Pi = cos^f[{ml-s){-4mj+ml-s)Co + 2s{Bo{s)-Bo{ml))+2{s-ml)], 

f/ = -sinffGo, (4) 


satisfying the gauge invariance ■ {f( g^'' — FIk'^q^) = F^k^ ■ [k■ qg>^^ — k''q>^) = 0 with k-q = (s —m ^)/2 
and e^^''°'^kfj_ki3qa = 0. In Eq. Q, Gq = Go{0,s,ml,m‘j,m‘j,m‘j),Bo{s) = Bo^s^m^mj), Bo{ml) = 
i?o(m^, TO^, TOj) with the analytical expressions being given in Appendix [A| 

The IT-loop contribution functions Ff’’^ and Ff'^{i = 1,3) in the Feynman gauge can be found 
in [in [50] and are collected in Appendix while Ff)’^ = F^’^ = 0. As first pointed out in m the 
lE-loop contribution to the hZ^ vertex diagram contains divergence which is cancelled by the divergent 
part of the IT-loop contributions to the Z — 7 mixing, while fermion-loop contributions to the Z — 7 
mixing are proportional to the final photon momentum squared and therefore vanishes for photon being 
on-shell. 

Following with the notations in Ref. [53], the amplitude Ai can be decomposed as 


M = 


1 

IGtt^ 


E E 

2=1,2,3,6 a=± 


A“Gf 


(5) 


and the form factor Gf" sum all the diagrams 


G± = G7± + Gf ± + G7± + g7°"± , ( 6 ) 

where the form factors G^^, Gf ^, G®^, G^™^ denote the contributions of the /177 vertex diagrams, the 
/ 1 Z 7 vertex diagrams, the t-channel hee vertex corrections and the box diagrams, respectively. Contri¬ 
butions from Z — "f mixing diagrams are also included in Gf ^. The matrix elements are given by 

Af = v{p+)i^Pl5)li'ig^''k ■ P- - k''p^)u{p-)elik), 


^ The Z/electron box diagram is obtained by replacing W boson and neutrino with Z boson and electron, and attaching 
the photon to electron line since there is no 'yZZ interaction at tree-level. 
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FIG. 1: Representative diagrams for e+e“ —>■ h'y. (o)( 6 ) denote the s-channel vertex diagrams, (c)(d) 
denote t-channel vertex diagrams and Wj neutrino box diagrams, respectively. The CP-violating htt 

coupling is marked with the blue dot. 


Af = u(p+)(l±75)7i/ff''''M(p-)e*(fc), 

Ag = i;(p+)(l±75)7i.e^''"'^9a^/3M(P-)e*(fc), (7) 

where p± and k are the four-momenta of and the final photon, e* (fc) denotes the polarization vector 

of the final photon and the Levi-Civita tensor rpj^^ correspondences between and 

are 


/^ 7 ± _ /^ 7 ± _ 77/^7 
^1 - ^2 - 


■y-yzt 


e 

2s' 

e 


^3^ = -^AGj + '-^Gj), 


/~y7± _ ^ /^7 

^6 ~ ^^ 6 ’ 




z± 


2 s 

= Cf ± = - 


ez" 


= 


riZ± 

^6 


4swciy(s — m| -b imz^z) 

± a . 

(Gf+ - 


G 


3 J 


ez 


4swCiy(s — m| + imzT z) 
ez^ 

4swciv(s — m| -b imzTz) 


'-gD, 


G, 


6 > 


( 8 ) 


where s = z"*' = py + g’^, z~ = Py ~ Pa with = — 5 + 2 sin 9w, Pa = ~ 2 ’ 

mass and width of Z boson, respectively. The t-channel and box diagram contributions are displayed 
in Appendix]^ The gauge invariance of the e^e“ -A hy amplitude is maintained after summing over 
IF-loop contributions to Gg^, G^^ ar id Gg ^, Gg™^ [53] (the fermion-loop contributions to Gg^, G^^ 
are zero), namely the relations in Eq. (BI 6 I. 

The unpolarized differential cross section of the process e“'"e“ -A h'y is 


da 


s-mt 


dcosO 647rs (IbTr^)^ 


K(|G+|^ + IGfl") + t\\G+\^ + \G^\^)+CVV} 


(9) 


where G^ are defined in Eq. (| 6 |) and 


CVV = - 


{s-mlf 


[(J-l_-J-l + )(l + COS0)2-(J-2_-J-2+)(l-COS0)2-(|Gg-p + |Gg+n(l + COs2 0)] (10) 


2 













5 


with 


■J'l± — 2 j '^6 ^1 ^6 ^1 i; — 2 j ''^6 ^2 *-^6 ^2 h 


( 11 ) 


The kinematic variables s, t and u are defined as s = {p+ + p-^, u = (p_ — A:)^ = — ^(s — m^)(l + cos0), 
t = (p+ — /c)^ = — 1 (s — TO^)(1 — cos 9) and 9 is the scattering angle between the directions of final photon 
and initial positron in the c.m. frame. The combination of the first two terms in Eq. is symmetric 
in cos 6* even in the presence of CP-violating hit coupling. While the last term, i.e. CW, can induce a 
forward-backward asymmetry which can be expressed as 


Afb 


ap — ctb 
CTp + as 


( 12 ) 


and 


ap 



da 

dcos9 


dcos9, Up = 



da 

dcos9 


d cos 9, 


(13) 


From the expressions in Eqs. § @ ([n]), the s-channel vertex diagram contribution to CW is like 
a + b cos 9 + c cos^ 9, where a,b,c are generic coefficients. Furthermore, b is proportional to the difference 
of right- and left-handed couplings, i.e. the axial vector coupling g’^ of the Z boson to e+e“. For 
nonzero tFi± and tF 2 ±, both the presence of CP phase ^ and strong phase are required, the latter of 
which arises from the threshold effects with on-shell intermediate W boson, Z boson (Z/electron box 
diagram) and/or fermions in the loops. In Fig. we show the energy dependences of loop functions of 
the s-channel diagrams. It is apparent that for 2mw < y/s < 2mt, the imaginary part may arise from 



FIG. 2: The energy dependences of various loop functions. For illustration, we have multiplied and 
Fg with a factor of 4:Q^Ncmf/m'^ ~ 25 and chosen ^ = 0,7r/2, respectively. 


Im(F’ 3 ’^) and Im(F^’'^), while for yfs > 2mt, ImFg and ImF/ contribute as well. From Eqs. ^ (|^, 
the form factors of s-channel vertex diagrams 




^z± e 
A'l = 


2s sw 


n,w 


2 

+ ‘2QtNc — ^{—2Qt + 


2s^c^(l 'm-'zis) 
■,± 


pZ,W) 


2s^c^(l m'^ls) 




(14) 


where the width of Z boson and light fermion contributions are not included here. In Fig. we show 
the real and imaginary parts of (7/ -I- Cf ^ for different There is a dip at i/s = 350 GeV for the real 
part distributions (a) (c) in the SM due to the destructive interference between the W boson and t quark 
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\^[GeV] ^[GeV] 

FIG. 3: The real and imaginary parts of for different 


contributions. As ^ increases, the values of real parts become larger since the cancellation is spoiled. For 
^ = 7r/2, only W boson loops contribute. From (b)(d), the imaginary parts are identical for different ^ 
as ^/s < 2mt since only W boson in the loop is on-shell, and differ as ^/s exceeds the tt threshold: for ^ 
varying from 0 to 37r/5, the imaginary parts increase. 

The box diagrams have all partial waves, i.e. cos" d,n > 0, in general. The form factors and 

C 2 °^^ of box diagrams depend on the scattering angle 0 and can provide an imaginary part if the W, Z 
bosons are on-shell. From Eq. Q and Fig. 11 the s-channel vertex diagram contribution interferes with 
the box diagram contribution destructively in the SM. In the presence of CP-violating hit coupling, the 
total cross section and differential cross sections for various ^ are discussed in section |III| 

The relation between Apg and CP violation in this process can be shown via the symmetry of the 
helicity amplitude A4\^r, where A/2 = A_ = —A+ with A_ and A+ being the helicities of initial electron 
and positron respectively, and r = ±1 are the helicities of the final photon. Since we have neglected 
the electron mass, conservation of the the electron chirality leads that the helicities of initial electron 
and positron are opposite. Under parity (P) transformation all helicities change signs, while charge 
conjugation (C) switches A+ and A_ (exchanges a particle with its antiparticle), that is A —)■ —A. Thus 
P, C and CP invariances give rise to |44l ^ 


Mx,r{0)=M^X,-r{9) 

Mx,r(0)=M-x,r(7r-0) (15) 

Mx,r(0)=Mx,-r(7r-0) 


respectively, up to phases which are however unimportant here |18j . For unpolarized beams and the final 
photon polarizations being summed. 


\Mx,t\'^ = |Al++(cos0)p -I- |Ad+_(cos0)p -I- |Al_+(cos0)p -I- \M _(cos0)p. (16) 

A,r=±l 


^ We would like to thank the authors of Ref. m for conversation. 
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FIG. 4: Left: cross sections as a function of c.m. energy ^/s for the CP phases ^ = 0 (black, solid), 7r/4 
(red, dashed), tt/2 (blue, dotdashed), Stt/S (orange, dotted). Right: cross sections as a function of cos^ 
for the c.m. energies ^/s = 250 GeV (blue, solid), 350 GeV (red, dashed), 500 GeV (green, dotdashed). 


(17) 


P, C and CP invariances imply the relations 

|AIa,tP = 2(|AI++(cos6»)p + |AI+_(cos6')p), 

A.t=±1 

|AIa,tP = |AI++(cos0)p + |AI+_(cos0)p + |AI++(-cos0)p + |AI+_(-cos6>)p, (18) 

A.r=±l 

y] |MA,rP = |M++(cos0)|^ + |M_+(-cos0)|2 + |M++(-cos0)|2 + |M_+(cos0)|2, (19) 


A.t=±1 


respectively. From Eqs. (18) (19), if C and CP are conserved the matrix elements squared are symmetric 
in cos 9. Therefore, a nonzero forward-backward asymmetry indicates both C and CP violation. 


III. NUMERICAL RESULTS 


In this section, we will give the numerical results of total cross section and differential cross sections 
at typical c.m. energies of future e+e” colliders [HI]. In practice we insert the Feynman rules of Eq. ([^ 
into FeynArts-3.9 [IS] model file, and calculate the amplitude automatically using FormCalc-8.4 [30] and 
LoopTools-2.12 [3^. We have checked both analytically and numerically the ultraviolet (UV) finiteness 
of the amplitudes. 

Fig. 0 shows the cross sections of e“''e —>■ h'y for different c.m. energies ^/s and CP phases From 

the left plot, we can see that the cross section grows as ^ increases from 0 to 37r/5. The shift of c.m. 
energy corresponding to the maximal value of the cross section from 250 GeV to 350 GeV is due to 
the competition of the W-loop functions with the t-loop function of the s-channel diagrams and their 
interference with the box diagram contribution. From the right plot, we can see that the cross section 
at a/s = 350 GeV drops rapidly (~ 10 times) as |^| decreases. Thus measuring the cross section of 
this process at 350 GeV can be helpful to search for CP violation in the hti interaction. Besides the 
total cross section, the differential cross section is also important. Fig. shows the normalized angular 
distributions for different a/s and The distribution is symmetric in the SM (^ = 0) and becomes 
asymmetric in the forward (cos 9 > 0) and backward (cos 0 < 0) regions as ^ deviates from zero, even in 
the case of f = 7r/2 Since the coefficients Vi± and J^ 2 ± in Eq. (11) depend on the scattering angle 9, the 


distributions for ^ yf 0 are not parabolic in shape. The asymmetry is most apparent at a/s = 500 GeV. 


^ This depends on our assumption of the couplings hWW and hZZ. In the 2HDMs [62], however, the CP-odd Higgs does 
not couple to the gauge bosons VP, Z at tree-level. 













FIG. 5: Normalized distributions of the scattering angle 9 for different ^ and y/s. 6 is defined as the 
polar angle between the momenta of final photon and positron. For ^ = 0, the distribution is 
symmetric. For ^ = 7r/2, ^ = 7r/4 , Stt/S, the distributions are asymmetric. 


Concentrate on the fist two terms in Eq. (101 and define the ’’differential cross sections 


dacpyj'^er) 

dcos9 

dacpvO^ox) 

dcos9 

dacpv 

dcosd 


= i4:{[(Cg-*C'™‘- - C'e-C™'-*) - (Cg+*C'™’'+ - C'g+C'™^+*)](1 + cos 6»)2 

- [(Cg-*C'™‘- - Cg-Cr*-*) - (Cg+*C'™’'+ - C'g+C™^+*)](1 - cos6»)2}, 

= K {[{ Cq * C ^°^- - - { C +* C \°^+ - C ' 6 + C ' J ^™+*)](1 + cos 6»)2 

- - C^Cr) - - C'g+C 2^™+*)](1 - cos0)2}, 

= i^{[(Cg-*Cr - Cg-Cr*) - (Cg+*C+ - Cg+C+*)](1 + cos0)2 

- [{CrC^ - C^Cr) - {C+*C+ - C+C^m - cos9)% 


( 20 ) 

( 21 ) 

( 22 ) 


with 0^2 = 1 ^ 1 ( 2 ^ + C\° 2 ^ and the overall factor K = i{s — m^)^/(4®7r®s). For 

the s-channel vertex diagrams thus ducpy (ver)/d cos 0 is proportional to cos0. For the 

box diagrams, the dependence on cos0 is not manifest since the expressions of g}™^ and G^™^ are 
complicated as displayed in Appendix]^ In Fig.[^ we show dcrcpv(vertex)/dcos0, dcrcpv(box)/dcos0 
and dacpv/dcos9 numerically for different CP phases ^ at y/s = 250 GeV, 350 GeV, 500 GeV. From the 
plots in the first two rows, the s-channel vertex diagram contribution and box diagram contribution to the 
difference of cross sections in the forward and backward regions, i.e. crp —crp are destructive. Furthermore, 
the magnitudes of both are maximal for ^ = 7r/2, since G^ are proportional to sin^. At y/s = 250 GeV, 
the magnitudes of dcrcpv (ver)/d cos 0 and dacpv O^ox)/d cos 6 are comparable so that there is a large 
cancellation between them and the resulting differential cross section dacpv/d cos 9 is small. Besides, 
cTp — cTp at y/s = 250 GeV is always positive. At y/s = 350 GeV, both dcrcpy (vertex)/d cos 0 and 
(^CTcpy (box)/d cos 0 are enhanced, since the loop function C/ gets its largest value near the tt threshold. 
The resulting distribution of dacpv/d cos 9 is not monotonic and the difference ap — ctb is positive for 
^ = 37r/5 while negative for ^ = -n 1/1:12 after integrating dacpv/d cos 9 in the forward and backward 
regions. The magnitude of the box diagram contribution at y/s = 500 GeV is about 3 times larger 
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|-^=250GeV| 


''A 

....f=;r/4 

-^=7t/2 

-^=37tI5 






-1.0 -0.5 0.0 0.5 1.0 

cos9 



FIG. 6: The distributions of dcrc’py(vertex)/dcos0, dcrcpy(box)/dcos0 and dacpv/dcosO at c.m. 

energies ^/s = 250 GeV, 250 GeV and 500 GeV, where we have not included the light fermion 
contributions. The dotted, dashed and solid curves correspond to ^ = 7r/4, •7r/2, 37r/5, respectively. 


than the box diagram contribution at yfs = 250 GeV while the s-channel vertex diagram contribution at 
^ = 500 GeV is smaller. As a result, dacpvId cos 9 ai ^Js = 500 GeV is two orders of magnitude larger 
and the difference ap — ap is negative. 

To illustrate the asymmetry quantitatively, we present the values of App in Fig. and Table |T] {App 
before cuts are applied) given the c.m. energy ^/s and GP phase The App depends on both the 
difference ap — ap and the total cross section, so App and up — up have the same sign. We find that 
the magnitude oi App tends to be small at low ^/s and get larger as y/s increases®. In Table {App 
before cuts are applied), the App is positive at y/s = 250 GeV but negative at y/s = 500 GeV, and the 
latter is two orders of magnitude larger since ap —ap at= 500 GeV is two orders of magnitude larger 
while their total cross sections are comparable, see Fig. Furthermore, the App at y/s = 350 GeV with 
^ = 7r/2,37r/5 are tiny in consideration of the small ap — ap and the enhancement of total cross sections. 


® For illustration, we also show App y/s = 400 GeV for varying cos^ in Fig.j^ 
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FIG. 7: Afb generated for various CP phases ^ and c.m. energies ^/s = 250 GeV (blue, triangle), 
350 GeV (red, triangle), 400 GeV (purple, disk), 500 GeV (green, rectangle). 
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IV. COLLIDER SIMULATION 


In this section, we will simulate the signal process e+e“ —>■ /17 —>■ 667 and its backgrounds at the future 
high luminosity {C = 1 ab“^, 3 ab“^, 10 ab“^) e“e+ colliders at the typical c.m. energy y/s = 250 GeV, 
350 GeV, 500 GeV. To generate signal events, we obtain the amplitude squared for e+e” —)■ /17 with the 
help of FeynArts [3S] and FormCalc/LoopTools packages [15] and then pass them to MadGraph [551151] . 
The background matrix element for e+e“ —>■ 667 is generated directly using MadGraph. The entire 
analysis is done at the parton level, with the following event selection cuts being applied: 

(а) pf > 20 GeV, pj > 25 GeV, < 2.5, /^Rbb,b-,,b^ > 0.4, (23) 

(б) \m^,l - m?,| < 15 GeV, x 0.5% (24) 

where plp and ry* denote the transverse momentum and pseudo-rapidity of the particle i, respectively. 
The spatial separation between the objects k and I is denoted by ARki- The photon in the signal event 
exhibits a harder transverse momentum to balance the momentum of Higgs boson than the photon in 
the background which is mainly radiated from the initial electron and positron and peaks in the small 
Pt owing to the collinear enhancement. The 6 -tagging efficiency Cb and the mis-tag probabilities Cj^b for 
light jets in our analysis are [55] 

eb = 0.9, ec^b = 0.1, ej^b — 0, for j = u, d, s, g, (25) 


and at least 1 6 -jet is tagged. 

Fig. 1^ shows the angular distributions of the signal processes e+e“ —>■ /17 —667 with the CP phases 
^ = 0, 7 r/ 4 , 7 r/ 2 , 37 r/ 5 . The p^ and cuts in Eq. (231 will constrain the maximum value of 0 and hence 
reduce the heights of the peak in the forward/back^rd region, while the other cuts in Eqs. (23) (24) 
are expected to have less impact on the distributions. The forward-backward asymmetries of the signal 
with different c.m. energies at a given value of ^ are displayed in Table [T] We find that the change after 
the cuts are applied is substantial at y/s = 350 GeV, including the signs. This is because the Afb 
at y/s = 350 GeV is very small and the sign is sensitive to the difference of events in the forward and 
backward regions. To illustrate the impact of the cuts on the Afb, h is helpful to write Afb = ANg/Ng 
as in Eq. (Cl), where ANg is the difference of the numbers of events in the forward and backward 
regions, and Ng is the total number of events. The angular distributions at y/s = 250 GeV are nearly 
symmetric in Fig. [^ so the cuts are expected to have more impact on Ng rather than on ANg, and Afb 
will become larger in magnitude after the cuts are applied. On the other hand, the angular distributions 
at y/s = 500 GeV are apparently asymmetric. So the cuts will have more impact on ANg and the 
magnitudes oi Afb are reduced. For y/s = 350 GeV, the distributions are moderately asymmetric, and 
ANg is very sensitive to the curve shapes, so that the Afb are greatly affected by the cuts. Fig. [^ shows 
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FIG. 8 : Angular distributions of the signal process = 0, 7 r/ 4 , 7 r/ 2 , Stt/S) that have passed the event 

selection cuts. 



COS0 


FIG. 9: Angular distributions of the SM background (bkg) that has passed the event selection cuts. 

The distribution is symmetric. 


the angular distributions of the SM background process e“'"e“ —>■ 667 , which remain nearly symmetric in 
cos 9 [ 66 ] with the cuts bein g applied. 

In Table [n] and Table III we show the outflows of the cross sections of the signals and background (bkg) 
at y/s = 350 GeV, 500 GeV and the significances corresponding to the integrated luminosity C — 1 ab“^, 
3 ab~^, 10 ab~^. If there exists GP violating hit interaction, the cross section at yfs — 350 GeV is 
strongly enhanced. On the other hand, the magnitude of Aps at \fs = 500 GeV is maximal which will 
be appropriate for the achievement of a larger statistical significance of the asymmetry. The foremost 
SM bkg is from the process e+e” —>■ —>■ 667 and the behavior of its cross section is ~ 1/s [Ml- From 

the tables, we see that the bkg cross section is larger at yfs = 350 GeV than at yfs = 500 GeV before the 
photon energy resolution is required. Fig. [T 0 |shows the signal and bkg distributions of the invariant mass 
of bi) and the photon’s recoil energy following the basic selection cuts (a) in Eq. (23). So the mass window 
cut ~ 15 GeV and the requirement of the photon energy resolution are efficient to improve the signal 
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WbblGeV] 



FIG. 10: Left: the invariant mass distributions of the b jets. Right: the distributions of the photon’s 
recoil energy. The signal is in red, while the background is in blue. The events are collected for the 
integrated luminosity £ = 10 ab“^ and the c.m. energy = 500 GeV. 


significance. We have chosen the planed resolution AE-y/E^ ^ 1.7% [SS] and the optimistic resolution 


/ E^ ~ 0.5% [SS] in our analysis, and the remaining cross sections are displayed in Table and 


AEy , _ 

Table III We can find that the photon energy resolution has more impact on the bkg at y/s = 350 GeV 
than at y/s = 500 GeV, due to the fact that E^ at y/s = 350 GeV is smaller. For the planned photon 
energy resolution, the signal significance Sj'/B can reach 3 (t with the integrated luminosity £ = 3 ab~^ 
and about 5cr with £ = 10 ab“^ at ^/s = 350 GeV, and only 3cr with £ = 10 ab“^ at ^/s = 500 GeV. For 
the optimistic photon energy resolution, the significance increases to more than 3 (t with £ = 1 ab~^ and 
5cr with C — 3 ab“^ at i/s = 350 GeV, and more than 2.1cr with £ = 3 ab“^ and 4cr with £ = 10 ab“^ 
at Vs = 500 GeV. 


TABLE I: The forward-backward asymmetries of the signal before and after the cuts in Eqs. (231 (24) 
for different ^ and ^/s. 


ApB before cuts applied 

Vs = 250 GeV 

Vs = 350 GeV 

Vs = 500GeV 

II 

0.0392 

-0.0231 

-0.5478 

C = V2 

0.0417 

-0.0018 

-0.4575 

^ = 37r/5 

0.0365 

0.0029 

-0.3987 

ApB after cuts applied 

Vs = 250 GeV 

Vs = 350 GeV 

Vs = 500GeV 

II 

0.0464 

0.0058 

-0.5090 

II 

0.0596 

0.0136 

-0.4312 

^ = 3n/5 

0.0460 

-0.0072 

-0.3616 


TABLE II: Cutflow of the cross sections (in fb) of the signals with ^ = 0,tt/A,tt/2,3tt/5 and the back- 


,-1 


ground at ^/s = 350 GeV with the integrated luminosity £ = 1 ab 3 ab ^ and 10 ab The SM value 
BR{h —>■ bb) = 57.7% and the b-tagging efficiency and the mis-tag probabilities in Eq. (25) have been 


included. 
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yi = 350 GeV 

II 

0 

signals (5) 
f = 7 r /4 ^ = ty/2 

^ = 37 r /5 

bkg (B) 

cuts (a) in Eq. ( 

23 

1 

0.017 

0.053 

0.121 

0.144 

477.55 

1771^5 — iTihl ^ 15 GeV 

0.017 

0.053 

0.121 

0.144 

14.70 

AE.y <E.,x 1.7% 

0.017 

0.053 

0.121 

0.144 

6.65 

AE.y <E.yX 0.5% 

0.017 

0.053 

0.121 

0.144 

1.93 

S/B 

0.003 

0.008 

0.018 

0.022 


Sj^fS with 1 ab~^ 

0.212 

0.645 

1.488 

1.762 

6.65 

Sjy/B with 3 ab~^ 

0.367 

1.118 

2.577 

3.052 


Sj^fB with 10 ab 

1 

0.669 

2.041 

4.705 

5.573 


SIB 

0.009 

0.027 

0.063 

0.074 


SjVB with 1 ab“^ 

0.393 

1.197 

2.760 

3.269 

1.93 

SjV~B with 3 ab~^ 

0.680 

2.073 

4.781 

5.662 


SjVB with 10 ab 

1 

1.241 

3.785 

8.728 

10.338 



TABLE III: Cutflow of the cross sections (in fb) of the signals with ^ = 0, 7 r/ 4 , 7 r/ 2 , Stt/S and the back- 


,-1 


ground at yfs = 500 GeV with the integrated luminosity £ = 1 ab~^, 3 ab~^ and 10 ab“^. The SM value 
BR{h —)■ bb) = 57.7% and the b-tagging efficiency and the mis-tag probabilities in Eq. (25) have been 
included. 


= 500 GeV 

II 

0 

signals {S) 

^ = 7 r /4 ^ = tt/2 

5 = 37 r /5 

bkg (B) 

cuts (a) in Eq. ( 

23 

1 

0.026 

0.041 

0.067 

0.073 

225.53 

\mjji — iTihl S 15 GeV 

0.026 

0.041 

0.067 

0.073 

6.12 

AE.y <E.yX 1.5% 

0.026 

0.041 

0.067 

0.073 

5.86 

AE.y <E^x 0.5% 

0.026 

0.041 

0.067 

0.073 

2.90 

SjB 

0.004 

0.007 

0.011 

0.012 


Sj^/B with 1 ab“^ 

0.337 

0.540 

0.875 

0.951 

5.86 

Sj^/B with 3 ab“^ 

0.584 

0.936 

1.515 

1.647 


SI\fB with 10 ab 

1 

1.067 

1.709 

2.766 

3.008 


SjB 

0.009 

0.014 

0.023 

0.025 


Sj'/B with 1 ab~^ 

0.480 

0.769 

1.244 

1.353 

2.90 

Sjy/B with 3 ab~^ 

0.831 

1.332 

2.155 

2.344 


Sjy/B with 10 ab 

1 

1.517 

2.431 

3.935 

4.279 



In order to estimate whether the ApB can be measured at the future high luminosity e+e“ colliders, 
we calculate the significance of the expected asymmetry with which a particular CP phase ^ would 
manifest nni 


S = \Afb\ 




y/ Ns + Nh 


(26) 


where Ng, N), are the n umb er of signal and background events, respectively, and ApB is the theoretical 
asymmetry given in Eq. ( 12 ). The significances at ^/s = 500 GeV are larger than those at ^/s = 350 GeV. 
For the c.m. energy ^/s = 500 GeV and the integrated luminosity £ = 10 ab~^, the significances are 
given by ^ = 1.23cr, 1.68cr, 1.57cr for ^ = 7r/4, 7r/2, 37r/5, respectively. 


V. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have investigated the effects of the CP-violating hit coupling in the process e+e“ —>■ 
£ 7 . Our numerical results show that the cross section can increase significantly for the allowed CP phase 
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^ and the c.m. energy. For -y/s = 350 GeV and ^ = 37r/5, the cross section is about 10 times of that in 
the SM. A preliminary simulation for signal e“'"e“ —>■ hy —>■ 667 and its backgrounds has been carried out. 
For y/s = 350 GeV, we can observe the signal at about 5tT with the integrated luminosity £ = 3 ab~^ 
with ^ G [7r/2, 37r/5]. For ^/s = 500 GeV, we can observe the signal at more than 2.1cr with £ = 3 ab~^ 
with ^ G [7r/2, 37r/5]. For £ = 10 ab“^, the significance can be greater than 5cr for ^/s = 350 GeV and 
than 4(7 for 500 GeV with ^ G [7r/2, 37r/5]. 

Compared with the Aps in the Higgs decay h —)■ the App can be greatly enhanced in the 

production process. App can reach -0.55 for ^ = 7r/4 and y/s = 500 GeV. Due to the large backgrounds, 
the significance of the expected App can be observed at only 1.68 ct with £ = 10 ab“^ and y/s = 500 GeV. 
We should emphasize that it is essential to trigger the single photon in the final state to separate the 
bottom jets arising from scalar or vector bosons, in order to isolate the signal from backgrounds, especially 
for high c.m. energy. 
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Appendix A: Analytical expression of scalar functions 


In this appendix, we give the analytical expressions of two-point and three-point scalar functions Bq 
and Cq. The scalar functions Bq and Co are defined as m 


Bo{pl,ml,ml) = £*^( 271 ^)'^ ^ f d^k 
Co{pl,pl,pl,rnl,rnl,rnl) = f d^k 

and the analytical expressions are [Ml |69] 


1 


— m\][{k -I- pxY — TO2] ’ 

1 


- m\][{k +piY - rri^[{k + pi + P2)‘^ - Wg] ’ 


(Al) 


Bq{ 


2 2 
s.m ,m 


) = A, - 

Jo P 


= A,. - log — -k 2 - 25 (A), 
k 


BQ[m\,mJ ,7X1^) = Aj — / dxlog 

Jo 


—x{\ — x)'m\ + xm\ + (1 — x)7 




= A,-log 2 - 25 ( 7 ), 


Co (0,0 


^2 _2 _2 _2 


, m ) = — 


) = - fdx f \y 
Jo Jo 


— xy — X y) -\- m? — is 




Co(0,0, s, m^, TO^) = — dx dy 

Jo Jo 


— xy — X y) m? — is 




ol — x 


Co(s,0,m^,m^,m^,m^) = - dx dy- - 2 I-^-TV-TV7- 2 -^ 

' ' (s — mf^jxymf^x{x — 1)— IS 


Jo Jo 
1 -2tA 
4m^ A — T 


(/(t)-/(A)) 


(A2) 
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where = | 
functions m 


7 b + log 47 r, e = 4 — £), 7 b = 0.5772 • • • is the Euler constant, r = A = and the 

nih s 


/(t) 


arcsin^ t > 1 , 

-iVl-'r(log^ -i7r)2 0 <r<l, 


(A3) 


with r]± = 1 zL -v/l — r. 




•\/'r — 1 arcsin^ -s/r 
JVl-T(log^ -*7r) 


T > 1, 

0 < T < 1, 


(A4) 


Appendix B: Loop functions and form factors 

In this appendix, we will give the detailed s-channel hE-boson loop functions and t-channel and box 
diagram form factors of the process e“'"e“ —>■ ^ 7 , computed in [S^- 


1. s-channel diagram contributions 

W boson loop functions 

F^’W = (Z(^+6)(4C24-Bi3)+m^(-7Co + Cn+Ci2) 

+ s{5Co + Cii — C 12 ) — Bi2 + B 23 , 

pz,w ^ 1 Z!^( 1 _ 2 c^)( 4 C 24 - 5 i 3 )+™^C^( 7 Co-Cii-Ci 2 - 2 C 23 ) 
2 m^ 

4" + Cii — C 23 ) + 2m^(c^ — l)Co 

~ c^[s(5C'o + Cii — C 12 + 2(721 ~ 2(723) + 32(724 — i?i 2 — 6 B 13 — 523 ] 

+ 2 5(1^0 ~ C*!! ~ 1^21 + <723) — Bi 2 + i ?237 

2 

F3’ = 2[—^(1 — 2c^)-l-2(1 — 6 c^)]((7i2 + ( 723 ) + 4(1 — 4c^)(7o, 

where the two-point and three-point functions are 

Bi2 = i?o(S7 TO^), i?i3 = i?o(w^, m^, TO^), 

B 23 = Bq(0, 171^,171^), (7o;ij = (7o;ij(s,0,TO^,m^,m^,TO^). 


(Bl) 

(B2) 


(B3) 

(B4) 


(B5) 


2. t-channel diagram contributions 


where 


-(7®=*= = 




Ct = 0, 


-I- crossed (i = 1, 2, 3), 

4c^ 


A^+ = (7i2(m2,f,m^,m^,0,m^), 

Af=^ = (z=^)^(7i2(m2,t,m^,TO|,m2,TO|), 
and the crossed terms are obtained by substituting t ^ u. 


= 0, Ar± = 

Af ± = Af ± = 0 


= 0 , 


(B 6 ) 
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3. box diagram contributions 


Cl 


box± 


[B'^ + crossed(_B2 )] + 


box± 


a 


_^boxib 


e'^TOvF r 

4s^ 


[B 


+ crossed(Bj^=^)] 


e "12 dZ± 
4,3 ^3 -“i > 

dZ± 

4„3 3 ^2 ) 

e^"»2 dZ± 

4,3 3 -°3 > 


where® 


C, 


VK+ 


box± 


= 0 , 


B. 


bY+ 


^w- 


= 4 (£)o + -Dll — Di3 + D23 — -D25), 

= 4(—-D 12 +-D 13 +-D 23 — D 26 ), 

= 2 + Di 3 + -D25 + D26 — D23 — D24) 

+ 7t(—-D 23 + -D 26 ) + m(3-Di3 + 7-D25 — 7 D 23 ) — 2 O-D 27 ], 

= 0 , 


R ^ 

-^ 1 , 2,3 

with the four-point functions 

Do-ij = Do;ij(mg,TOg,m^,0,s,u,m^,0,m^,TO^) 

and 


B 


— 2 (z^)^(—-Dll-I--D12-I-D22 — D24), 


±)2, 

= ( 2 ^)^[s(D 22 — -D 24 -I--D 25 — D 26 ) + 2 -D 27 ] 


i3f± = 2 (z±)2(D22-D26), 

jZ± 


with 


Do-ij = -Do;y(mg,TO^,mg,0,t,u,mg,TO|,TO|,mg). 

The crossed terms are obtained by substituting t ^ u. In fact, we find that 


= crossed(-Bf *), 


Z±\ 


= crossed(Cf“*), 


box± \ 


(B7) 

(B8) 

(B9) 

(BIO) 


(Bll) 

(B12) 


(B13) 

(B14) 

(B15) 


and the combination (7®™=*= 
In Fig. 


_ 1 _ (^boxi jg symmetric in cos 0. 
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we show the SM = 0) the differential cross sections da/dcos9 at c.m. energies y/s = 
250 GeV, 250 GeV and 500 GeV. The solid, dashed and dotted curves denote the full, box diagram 
and s-channel vertex diagram contributions, respectively. We can easily find that the box diagram 
contributions interfere destructively with the s-channel vertex diagram contributions and dominate at 
higher c.m. energy in the SM. 

Finally, we have checked explicitly the following relations using the LoopTools package |46) . 


cr if) = ct if) = ctif) = cl^if) = o> 

criw)+c^-m = ^, 

CT + = 0 , 

C'3'+(IF) + G3^+(IF) -b (71+ -b (7®“+ = 0, (B16) 

which are the results of gauge invariance, and the form factors Cf’^^{f) and G 3 ’'^+(IF) denote the 
fermion-loop and W-loop contributions to respectively. 


® There are several typos in the above expressions in Ref. m which have been corrected here. 
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FIG. 11: The differential cross sections da/dcosO at c.m. energies ^/s = 250 GeV, 250 GeV and 
500 GeV in the SM. The solid, dashed and dotted curves denote the full, box diagram and s-channel 

vertex diagram contributions, respectively. 


Appendix C: Significance of the forward-backward asymmetry 


In this appendix, we will give the detailed derivation of the significance of the expected asymmetry 
in Eq. (26), which has been discussed partly in in uni mini [201 HD .We define the theoretical asymmetry 
Afb and the measured asymmetry as 


^FB 

^meas 

/ipB 


ANs 
N s ’ 

AN 

Ns+Nb^ 


(Cl) 

(G2) 


where ANg = Nf — and Nf and Nf are the number of events in the forward and backward regions, 
respectively. Ng = Nf + N^ is the total number of signal events, and Nb is the number of the background 
events. AN ~ ANg if the background does not contribute to the asymmetry or the contribution is much 
smaller. Write Nf = F and N^ = B, we have 


The statistical error of ApB is expressed as the error propagation, 

2 fdNp'Q 2 2 I /^^FBn2 2 

For a Poison distribution, ap=F and = B. Thus we obtain 

1 _ /t2 
2 _ ^FB 

^ Ng ’ 


(G3) 


(G4) 


(G5) 


which is approximate cr^ = l/\/Ng for a small ApB- Similarly, for the measured asymmetry the 

statistical error is 


(^meas)2 


1 / /imeas\2 

^ ~ l^FB ) 

Ng + Nb 


1 

Ng+Nb' 


(G6) 


The significance of the expected asymmetry S is defined in uni 


or in Refs. [Hin] 



5 = 


AN 

^Ng + Nb ■ 


(G7) 


(G8) 
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